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Abstract 

This paper presents a measurement of J/ipjip' differential cross sections in pp colli- 
sions at = 1.8 TeV. The cross sections are measured above 4 GeV/c in the central 
region (\r]\ < 0.6) using the dimuon decay channel. The fraction of events from B de- 
cays is measured, and used to calculate b quark cross sections and direct J/tp,Tp' cross 
sections. The direct cross sections are found to be more than an order of magnitude 
above theoretical expectations. 



1 Introduction 

Charmonium production is currently the best way to study b quark production at CDF at the 
lowest possible transverse momentum of the b. While the signature and signal-to-noise for 
the J/ip are excellent, the actual conclusions regarding b production are strongly dependent 
on the fraction of the data sample due to b decays. One obvious way of determining this 



1 



fraction is to use the decay distance of the J/ip state. In contrast, the fraction of ip'^s from 
h decays is thought to be close to one[|l|]. Thus, in determining the spectrum of h quarks 
from a study of ip' events, one does not have to worry about the fraction of if)' production due 
to B decays. In the course of this analysis, however, a very large zero-lifetime component 
for the observed ip' signal has been seen. Therefore, one still has to make use of the lifetime 
information. 

2 Detector Description 

The CDF detector has been described in detail elsewhere 0. We describe here briefly the 
components relevant to this analysis. The CDF coordinate system defines the beam line to 
be the z direction. R is the radial distance from the beam line, and (f) is the azimuthal angle. 
A solenoidal magnet generating a 1.4 T magnetic field surrounds the two tracking chambers 
used. The Central Tracking Chamber (CTC) is a cylindrical drift chamber surrounding the 
beam line. 

The CTC contains 84 layers, which are grouped into nine superlayers, with 5 superlay- 
ers providing axial information and 4 providing stereo. It covers a pseudorapidity range 
of \'q\ < 1.1. The Silicon Vertex Detector (SVX) is a silicon microvertex detector. It con- 
sists of four layers 2.9 to 7.9 cm from the beam line and provides high resolution tracking 
information in the R — (j) plane. When combined with a CTC track, it provides an im- 
pact parameter resolution of (13 + 40/Pp)yum and a transverse momentum resolution of 

^(0.0009Pj^)^ + (0.0066)^. Primary vertices have a Gaussian distribution with a width of 
about 27 cm, while the SVX only reaches to ±25cm, so only about 60% of the CTC tracks 
can be augmented with SVX information. 

Outside the CTC are electromagnetic and hadronic calorimeters, which provide five ab- 
sorption lengths of material before the Central Muon Chambers (CMU). The CMU consists 
of four layers of limited streamer chambers. It is divided into 72 segments which cover 85% of 
the (p region for |?7| < 0.6. Muons with below ~ 1.4 GeV/c range out in the calorimeters. 

The CDF trigger consists of three levels. In the first level, both muons must be detected 
in the CMU. They must be seperated by at least 0.09 radians in and pass a slope cut in 
the R — (j) plane. The slope cut is effectively a cut on the transverse momentum. To pass the 
second level, a fast hardware tracker must find at least one track that points to the correct 
muon chamber. For the third level, the full CTC track reconstruction is done. Both tracks 
must be found, and extrapolate to within 4cr of the muon chamber tracks in R — z and R — cf). 
a is the expected multiple scattering error added in quadrature with the measurement error. 

3 Event Selection 

Further cuts are applied offline to produce a purer sample. The CTC track and the CMU 
segment are required to match better than 3cr in P — 0, and 3.5cr for the J/'?/' and 3(J for 
the '?/'' in R — z. Hadronic energy is required in the calorimetry tower that points to the 
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muon chambers in the ip' events. Both muons are required to have Pj^ > 2.0 GeV/c and 
one muon must have PjJ > 2.8 GeV/c. The dimuon is required to have Irj'^'^l < 0.6 and 
Pj^'^ > 4 GeV/c. Runs with known hardware problems were excluded. Because there are 
much better statistics in the J/ip, a tighter definition of bad runs was used there, resulting 
in a smaller integrated luminosity. 

For the J/ip, the CTC track is beam constrained. In the ip', the tracks are vertex 
constrained, using the SVX information if there are three or more hits. The of this 
fit is required to be less than 10, with one degree of freedom. The resulting invariant 
mass is used to define a signal region of 3.0441 < m^^ < 3.1443 GeV/c^ for the J/ip and 
3.636 < m^^ < 3.736 GeV/c^ for the ip'. Sideband regions of 2.9606 < m^^ < 3.0274 GeV/c^ 
or 3.1610 < m^^ < 3.2278 GeV/c^ for the J/ip and of 3.52 < m^^. < 3.62 GeV/c^ or 
3.75 < m^^ < 3.85 GeV/c^ for the ip' are also defined. 

Figure 1 show the mass distribution of the events passing these cuts. There are a total 
of 26533 ±175 J/ip events and 896 ± 94 ip' events. 



4 B Fraction 

For events where both muons have SVX tracks with at least three hits, we vertex constrain 
the tracks to measure Ixy, the projection of the decay length onto the J/ip,ip' transverse 
momentum. This is converted into the proper lifetime of the parent by cr = Ixy/HPT/m) ■ 
jpcorr-^^ where F™'^'" is a Monte Carlo correction factor that relates the J/ipiip' tranverse 
momentum to the transverse momentum of the parent. 

The background cr shape is measured from the sidebands. It contains non-Gaussian 
tails, so we parametrize it by 



Bix) 



[l-f^- /_) 1 exp + f . ^ exp f^] , a; < 



where /+ is the fraction in the positive exponential, A+ is the lifetime of the positive ex- 
ponential, /_ is the fraction in the negative exponential, A_ is the lifetime of the negative 
exponential, and a is the width of the central Gaussian. The signal region is the fit to the 
function ■ fsigicr), where 

fs^gicT) = f.^ck " B{ct) + (1 - h^ck) [(1 - Ib) ' R{ct) + fB-R^ E{ct)] 

where R{ct) = /^.-^i^ exp + (1 — fr)jexp is the resolution function and 

E{ct) = ^ exp (-^f^)- fback IS the background fraction, fs is the B fraction, ctq is the B 
lifetime, a is the width of the Gaussian, is the fraction of non-Gaussian tails in resolution 
function, and A is the lifetime of non-Gaussian tails. R -k E{ct) is the convolution of the 
resolution function with an exponential. We fix ctq to 438 fim, as found by the CDF inclusive 
B lifetime measurement. H. fback is calculated by 
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ij" Vjrfi V 1 C) 


jj ridctioii 


Above 4 


19.6 ± 1.5 


4-5 


13.3 ± 1.0 


5-6 


15.9 ± 1.2 


6-7 


21.0 ± 1.7 


7 — 8 


25.2 ± 2.1 


8-9 


25.2 ±2.2 


9-10 


26.9 ±2.5 


10 - 11 


34.3 ±3.3 


11 - 12 


31.0 ±3.5 


12 - 13 


32.4 ±4.1 


13 - 14 


42.1 ±5.7 


14- 15 


27.6 ±5.0 



Table 1: Differential B fraction and cross section (from all sources) of J/ip. Errors are 
statistical plus systematic. 



In the ip' the data is fit using one unbinned log-likelihood fit. fhack is not fixed, but a 
Poisson term is included for the likelihood that they fitted value fluctuates to the calculated 
value, a is calculated on an event by event basis. Using the event-by-event error introduces 
non-Gaussian tails into the total resolution function, so fr is fixed at 1. The systematic error 
is estimated in part by allowing fr and A to be fitted. In the J /ip, the signal and sidebands 
are fit in seperate binned fits, fhack is fixed and a, fr and A are fit. Figure 2 show the result 
of the fits. 

4.1 Prp dependence 

It is expected that the B fraction in the J/ipiip' samples rises with Prp. To measure this, we 
divide the tp' sample into three bins from 4 — 6 GeV, 6 — 9 GeV, and 9 — 20 GeV. The 
size of the last bin is neccesitated by the low statistics at high P^. We then repeat the above 
fitting procedure in each of the P^ regions. 

In the J /tp sample, the "spectrum method" is used to calculate the P^ dependence. This 
is done in three steps. First, the J /ip P^ spectrum is obtained. Then, the spectrum of events 
with CT > 250;um is found. The two are normalized so that the ratio of the total area is 
the inclusive B fraction found above. The B fraction is then obtained by dividing the two 
distributions. The results are shown in Tables |I] and |^. 
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Pt (GeV/c) 


B Fraction (%) 


Above 4 


22.8 ±3.8 


4-6 


12.3 ±4.9 


6-9 


29.5 ±7.0 


9-20 


39.3 ±7.9 



Table 2: Differential B fraction and cross section (from all sources) of ■?/''. Errors are statistical 
plus systematic. 



5 Acceptance and Efficiencies 

The J/ip,ip' differential cross section is defined by 



dP-r a- e- J Cdt- APt 

where N is number of J/ip-,ip' candidates in the bin, a is the detector and kinematic 
acceptance, e is the trigger and cut efficiency, / Cdt is the integrated lumonosity, and AP^ is 
the size of the bin. The acceptance was determined from several Monte Carlo data sets. 
The Tp' acceptance used ip"s generated flat in Prp and r/. The J/ip acceptance used several 
sets where h quarks were generated in different Prp ranges, in order to provide sufficient 
statistics in all regions. The h quarks were generated using the next-to-leading order QCD 
calculation]^ and MRSDO structure functions. The quarks were fragmented to B mesons 
using Peterson fragmentation with eh = 0.006 ± 0.002. A fast detector simulation was 
used on the events, and the kinematic cuts were then applied to the events. The J/ip,ip' 
acceptance is the ratio of events that survived this process to the number of generated events. 

The b — >• J/ip,i'' acceptance used events where b quarks were generated and forced to 
decay to J/ip^ip'. The J/iJj,ijj' was forced to have the momentum spectrum measured by 
The Prp of the b quark is described by P™". P^™ is defined as the Pp such that 
about 90% of the b quarks in our sample have Pp > P™". The b quark acceptance is then 
calculated by 

_ Arfc(p//^'^' > 4.0 GeV/c, Ir^-^/V^-V^'l < 0.6, |/| < 1.0) 
~ Nb{P^ > P™", \y''\ < 1.0) 

The lack of a Pp cut on the b in the numerator corrects for the fact that we are unable to 
remove those events with Pp < P^*". 

The efficiencies of the three triggers where studied with unbiased J/ip events, as were 
the muon quality cuts. The efficiency of the CTC track reconstruction is 98.9 ± 1.0%. The 
efficiency of the muon reconstruction is 98 ± 1%. 
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6 Systematics 



Systematic errors can arise in three distinct areas: measuring the Jli)ii)' cross sections, 
determining what portion of that cross section is from B decays, and extrapolating those 
results to a 6 quark cross section. 

6.1 J/'4^')'4^' Cross Section Systematics 

Wc calculate the uncertainty due to the trigger shape by varying the shape of the level 1 
and level 2 trigger efficiency parametrizations. This is an 8% effect in the J/ip and a 6.6% 
effect in the ip' . We estimate the effect of the acceptance correction in the ip' by varying the 
acceptance by the statistical error and by fitting the acceptance to a fifth degree polynomial, 
which is a 3.1% effect. In the J/'^, fits to different order polynomials were used, giving a 1.5% 
error. A 5% error is assigned due to the J/iIj,iIj' polarization based on a ~ 10% difference 
between the acceptance of maximallly polarized J/ip^s with extreme values of a = ±1. The 
1% error measuring the cut efficiencies is neglibable compared to other errors. There is a 
2% error associated with each of the level 3 trigger efficiency, the CTC tracking efficiency, 
and the muon reconstruction efficiency. We assign a 4% systematic to the uncertainty in the 
luminosity. 

6.2 B Fraction Systematics 

Potential sources of systematic errors arise from the B hfetime, the spread of the acr distri- 
bution, and the fitting technique. 

A single value for the cr resolution is assumed in the J/ip fit, while a wide range exists. 
The effect of this is estimated by applying a cut of acr < 60/im, which removes the long tail 
in the distribution, acr is related to the opening angle between the muons, which decreases 
at higher Pj,. Since the J/ip from B decays have a harder Pj. spectrum, this may bias the B 
fraction, which is lower with this cut. The 7% difference is used as the systematic error. In 
the i/)', the unbinned fit uses the event-by-event error, so this effect is missing. Varying the 
B lifetime by one sigma changes the B fraction by 0.002 or less. This is negligible compared 
to other contributions. 

Changing the fitting technique and the form of the fitted function in the ip' can change 
the B fraction by 7%, which is used as the systematic error. 

6.3 B Cross Section Systematics 

We assign a 1.3% systematic error due to the uncertainty in the rapidity spectrum in the 
range 0.6 < \y^\ < 1.0. We take an uncertainty of 4.6% due to the uncertainty in the NDE 
spectrum with different values of fi and A. We assign 5% based on half the difference in the 
acceptance calculated with Peterson epsilon set to 0.004 and 0.008. A systematic of 5% is 
assigned based on the uncertainty in the CLEO momentum distribution. We assign 

a 1% error because of the statistics in the b to ■0' Monte Carlo. 
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In the V'', the effect of the large bin sizes for fs is estimated by fitting the B fraction to 
a straight line, and using a bin-by-bin correction factor. 



7 Cross Section 

7.1 J/tJj,tJj' Cross Section 

The cross sections are shown in Figure 3. The integrated cross sections are 

a{pp J/i), p//^ < 4 GeV/c, \r]\ < 0.6) ■ Br{J/i) ^T) = 
29.10 ±0.19(stat) tHUsyst) nb 

a{pp ip', pf < 4 GeV/c, |r/| < 0.6) ■ Br{i)' /i+zi") = 
0.721 ± 0.058(stat) ± 0.072{syst) nb 

The cross section from B decays is extracted by multiplying the differential B fraction 
by the cross section. We can also obtain a prompt cross section by multiplying the cross 
section by one minus the B fraction. The results are shown in Figures 3 and 4. The prompt 
theory curves are from |11|. 



7.2 Inclusive b Cross Section 

We now convert these cross sections into an inclusive b cross section. We do so using the 
standard formula, 

a{pp bX, \r]''\ < 1.0, > Pp*") = 

/ 



1 1 1 



2JCdt- ab Br{b J/ip, ip'X) ■ Br{J/ip, ip' /i/i) 



N{Pt) 



a(P^) ■ e(P^) 



The factor of | results from the fact that the b or the b can decay into the J/ip^ilj' . We use 
branching ratios of Br{b ^ J/i^X) = 1.16 ± 0.09% 0, Br{b ip'X) = 0.30 ± 0.06%[|, 
Br{J/i; fifi) = 6.27± 0.20%[§ and Pr(^' /i/i) = 0.88 ± 0.13%|T§. Our resuhs are 

a{pp bX, \y^\ < 1.0, P^ > 6.0 GeV/c) = 12.16 ± 2.07 /ib J/ijj 
a{pp bX, \y''\ < 1.0, P^ > 7.3 GeV/c) = 8.24 ± 1.34 /xb J/^ 
a{pp bX, \y^\ < 1.0, P^ > 8.7 GeV/c) = 5.20 ± 0.83 /xb J/V- 
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a{pp bX, \y''\ < 1.0, > 5.9 GeV/c) = 6.12 d= 2.04 /xb ^' 
a{pp bX, l/l < 1.0, > 8.3 GeV/c) = 3.85 ± 1.23 /xb ^/;' 
The error is statistical and systematic errors, added in quadrature. 
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Figure 1: Mass distribution (J/V'j V'') 
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Figure 2: J jijj and Lifetime Distributions. In the ijj' plot, the dark region is the background 
shape. The slashed region is the B component, plus the background shape. The clear region 
shows a clear excess of prompt events. 
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Figure 3: Differential Cross Sections. The B (prompt) values have been artificially offset 
by 200(-200) MeV/c for clarity. 
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Figure 4: J/-^ Differential Cross Sections 
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Figure 5: Ratio of data to theory for non-B Cross Section 



